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Abstract. The qualitative and quantitative understanding of near-horizon gravitational dynamics 
in the strong-field regime represents a challenge both at a fundamental level and in astrophysical 
applications. Recent advances in numerical relativity and in the geometric characterization of 
black hole horizons open new conceptual and technical avenues into the problem. We discuss 
here a research methodology in which spacetime dynamics is probed through the cross-correlation 
of geometric quantities constructed on the black hole horizon and on null infinity. These two 
hypersurfaces respond to evolving gravitational fields in the bulk, providing canonical "test screens" 
in a "scattering"-like perspective onto spacetime dynamics. More specifically, we adopt a 3+1 
Initial Value Problem approach to the construction of generic spacetimes and discuss the role and 
properties of dynamical trapping horizons as canonical inner "screens" in this context. We apply 
these ideas and techniques to the study of the recoil dynamics in post-merger binary black holes, an 
important issue in supermassive galactic black hole mergers. 
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A CROSS-CORRELATION APPROACH: MOTIVATIONS AND 

OBJECTIVE 

The general problem here discussed is the qualitative and quantitative understanding 
of near-horizon gravitational dynamics in the strong-field regime of black hole (BH) 
spacetimes. This represents a challenge both at a fundamental level and in astrophysical 
applications. The setting of our discussion is that of classical spacetimes in General 
Relativity (GR), with a focus on astrophysically motivated problems in which we adopt 
a numerical relativity methodology. 

A natural strategy to the study of spacetime dynamics consists in extending to the 
general relativistic setting the Newtonian description of the dynamics of gravitationally 
interacting bodies, namely celestial mechanics. This has proved extremely successful 
in unveiling the physics of compact objects. However, such an approach also meets 
fundamental obstacles in the general dynamical regime of a gravitational theory in 
which i) a priori rigid structures providing canonical references (such as symmetries 
or preferred backgrounds) are genetically absent, and ii) where global aspects play a 
crucial role. Here we rather adopt a complementary approach in the spirit of a coarse- 
grained description of the dynamics in which we renounce to the detailed tracking of 
the geometry (trajectories) of given compact regions (objects), and rather emphasize the 
global/quasi-local properties of the relevant dynamical fields. In particular, we aim at 



capturing the functional structure of the latter through appropriate correlation functions . 

From a physical perspective, we focus on what is sometimes called the "establish- 
ment's picture of gravitational collapse" [1]. Our current understanding of this problem 
can be summarized in the following chain of theorems and conjectures: 

i) Singularity theorems [2, 3, 4, 5]: if sufficient energy is placed in a sufficiently 
compact spacetime region, light locally converges in any emitted direction, trapped 
surfaces form and a spacetime singularity develops in their causal future. 

ii) Weak cosmic censorship conjecture [6]: to preserve predictability, the singularity is 
hidden from a distant observer behind an event horizon giving rise to a BH region. 

iii) BH spacetime stability (conjecture): GR dynamics drive the system to stationarity. 

iv) BH uniqueness (theorem) [7]: the final state is a subextremal Kerr BH spacetime. 

The establishment's picture of gravitational collapse is an intrinsically dynamical pic- 
ture. Accordingly, we adopt a methodology able to cope with the construction and anal- 
ysis of generic spacetimes, namely in an Initial Value Problem approach. This provides 
a systematic avenue to the study of the qualitative and quantitative aspects of generic 
dynamical spacetimes, respectively addressed through the use of tools in Partial Differ- 
ential Equation theory and by the numerical construction of spacetimes. 

Our focus is on the point iii) above, in an attempt to gain insight into the specific 
manner in which GR drives the system eventually to stationarity. Remarkably, tools 
devised in the study of i) and ii) (namely, the characterizations of BH horizons) prove to 
be very useful to reach this goal. In this sense, a crucial outcome of the systematic 
numerical exploration of fully dynamical vacuum spacetimes in recent years is that 
the a posteriori description of gravitational dynamics is rather simple. This specific 
observation is the main point supporting the applicability of a coarse-grained approach 
to the analysis of generic spacetime dynamics. 



A cross-correlation approach to BH spacetime dynamics 

In order to make concrete the previous considerations, we formulate the following 
cross-correlation methodology [8, 9, 10], whose specific goal is the development of 
qualitative insights into spacetime dynamics in the strong-field regime, identifying the 
key elements leading ultimately to appropriate quantitative effective descriptions: 

i) Spacetime dynamics is probed through the cross-correlation of geometric quantities 
h mn an /z out defined at inner and outer hypersurfaces, respectively, ^ mVL and Jf? out . 

ii) Hypersurfaces ffl mVL and J^ u t are taken as test screens responding to bulk dynam- 
ics. Spacetime is then explored in the spirit of an "inverse scattering approach". 

In our near-horizon and asymptotically flat context, the BH event horizon $ and future 
null infinity provide natural choices, respectively, for 3ff[ nn and J^ ou i (cf. Fig. 1). 
Quantities h mxi {v) and h ou t(u) are then constructed on <§ and as functions, respec- 
tively, of appropriate advanced and retarded times. We note that the cross-correlations 
(as time-series) between h-^iy) and h out (u) require a gauge mapping between v and u. 








FIGURE 1. Sketch of the conformal diagram for a spherically symmetric gravitational collapse, illus- 
trating the scenario for the cross-correlation approach to the analysis of bulk spacetime dynamics. 



BH HORIZONS: GLOBAL VS. QUASI-LOCAL APPROACHES 



BH horizons play a crucial role in our discussion, offering in particular a model for 
the inner hypersurface ^i nn in the cross-correlation scheme. Two approaches to the BH 
notion can be considered, both contained in the the standard picture of gravitational 
collapse. The first one is guaranteed by weak cosmic censorship: the black hole region 
g$ is a region of no-escape not extending up to infinity. In a strongly asymptotically 
predictable spacetime with J~(j r+ ) denoting the causal past of we have 
SB = M — 7 _ (J zr+ ). The boundary of the black region is given by the event horizon 
$ = dJ H j& . This traditional characterization of BHs for asymptotically flat 

spacetimes involves global spacetime concepts. In particular, event horizons $ are 
teleological objects whose location requires the full knowledge of the spacetime in the 
future and that can develop in flat regions. Event horizons are therefore not adapted for 
probing the BH spacetime during its construction in an Initial Value Problem approach. 

Alternatively, instead of characterizing the BH as the region that cannot send signals 
to distant observers, we can approach it as the region where all emitted light rays 
"locally converge". This is made precise by the notion of trapped surface, crucial in 
the singularity theorems appearing in the standard gravitational collapse picture. Given a 
closed surface C M with area element dA = y/q d 2 x, its normal plane can be spanned 
in terms of an outgoing null vector i a and an ingoing null vector k a . The outgoing 
expansion 6^> measures the rate of change of dA in the lightfronts emitted from 
along the outgoing null direction l a . The ingoing expansion 0w is defined analogously 
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Then 5? is a (future) trapped surface [2] if < and 0w < 0. The limiting case 
in which one of the expansion vanishes defines a marginally trapped surface: 0^ = 0, 



9^ < 0. If there exists a notion of outer region, e.g. associated with an asymptotically 
flat region, then (future) outer trapped surfaces [5] can be introduced as 6^ < 0, with- 
out any requirement on 0^ k \ A marginally outer trapped surface satisfies then 9^ = 0. 
In this context the BH can be characterized in terms of the trapped region ST , namely 
the set of points in spacetime belonging to some trapped surface. In the present con- 
text of a cross-correlation analysis, rather than in the BH or trapped region themselves 
we are primarily interested in its boundary. The latter provides a hypersurface to be 
employed as an inner screen ^fnn- hi this sense, the so-called trapping boundary [11], 
namely the boundary of the trapped region 3 ', would be a good candidate for J^ nn . 
Unfortunately, we lack an operational characterization of such trapping boundary. This 
is in contrast with the notion of apparent horizon, namely the boundary of the trapped 
region contained in a given 3-slice E, characterized as a marginally outer trapped sur- 
face 9^ = 0. Motivated by these difficulties, trapping horizons (namely worldtubes of 
marginally outer trapped surfaces) were introduced [11] as quasi-local models for the 
BH horizon, in particular in an attempt to gain insight about the trapping boundary. 



QUASI-LOCAL HORIZONS AS INNER TEST SCREENS 



We review the properties of dynamical trapping horizons, stressing those aspects of 
special relevance in our discussion. Let us consider a closed orientable 2-surface 
embedded in a 4-dimensional spacetime {■J%,g a b)i with Levi-Civita connection V a . 
Regarding its intrinsic geometry, we denote the induced metric as q ab , with Levi-Civita 
connection D a , Ricci scalar 2 R and area form e ab (area measure dA). Again, normal 
outgoing and ingoing null normals are denoted as i a and k a , normalized as £ a k a = — 1. 
This leaves a (boost) rescaling freedom £' a = f£ a , k la = f~ l k a , with / a function on 5?. 
We need the following elements of the extrinsic geometry. The outgoing expansion 

given in Eq. (1), and the outgoing shear are expressed as 



9^= q - b v a e b , o<2 = <f a f b v c e d -W<> 
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(2) 



whereas a normal fundamental 1-form Cl a associated with i a is given by 

d£ ] = -k c q d c y d £ c , 



(3) 



provides a connection on T^y. The transformation of these quantities under a null 
rescaling are: 0(0 = f9 (e \ ojp = foj£ and Q.P = ft? +D a (\nf). Finally, we need 
to control the variations of the outgoing expansion 9^ along normal vectors v" G I 1 ./ 
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where k-M = -v a k b VJ b , a and j8 are functions on y and 8 V is the variation operator 
associated with a change in the underlying surface y (cf. [12, 13, 14]). 

A trapping horizon [11] is (the closure of) a hypersurface Jf? foliated by closed 
marginally outer trapped surfaces: Jt? = \J te ^y t , with 0^ = 0. The properties of Jif 
as a horizon are characterized by the signs of 0^ and 8 k 9^\ i) the sign of 0^ controls 
if the singularity occurs in the future (6^ < 0) or in the past (0^ > 0), and ii) the sign 
of 8 k 0W controls the (local) outer (8 k 9^ < 0) or inner (40 W >0) character of 
with respect to the trapped region. For BHs the singularity occurs in the future and the 
horizon is an outer boundary. Therefore quasi-local BH horizons are modeled by future 
outer trapping horizons (FOTH): 9® = 0, ©(*) < and 5*0^ < 0. 

We can define an evolution vector h a on Jf? characterized by: i) h a is tangent to Jf? 
and orthogonal to y t , ii) h a transports y t to y t+ § t : = 1> an d iii) h a is written as 
h a = i a - Ck a . We also define a dual vector x a = i a + Ck a orthogonal to Jf. The sign 
of C fixes the point-like metric type of JF: C > spacelike, C = null, C < timelike. 
Considering first the spherically symmetric case (C constant on y t ), the trapping horizon 
conditions 0^ = 0, 5/,0^ = imply, using Eq. (4) and Einstein equations 

a®o® ab + T ab Pi b 

8 k 9 W ^ ' Pj 

for an outer Jf? (i.e. 8 k 9^ < 0) under the null energy condition. Therefore FOTHs Jt? 
can be either null or spacelike hypersurfaces. The first case corresponds to the stationary 
regime (with its the isolated horizon hierarchy [15, 16]), whereas the second case leads 
to dynamical horizons (DHs) [17, 18]. Beyond spherical symmetry, one can pose the 
question if a MOTS section y t of a FOTH can be partially spacelike and partially null, 
i.e. if it can happen C > in a part of y t and C = in another part. The answer is in the 
negative: transitions from stationarity to the dynamical regime happens "all at once" in 
sections y t of Jtf. This follows from the trapping horizon condition written as [use (4)] 

5h Q(e) = - 2 D C 2 D C C + 2^ C2 D C C-C8 k 9^ + 8 e 9^ =0 . (6) 

Under an outer condition 8 k 9^ < 0, a maximum principle can be applied to this elliptic 
equation to conclude that either C > (if 8 e 9 ^ = ojfj o ^ ° + T ab £ a £ b ^ at some point 

of y t ), or C = if and only if 8(9^ = everywhere on y t . In other words, it suffices 
that some energy crosses the horizon at a single point, for making the whole horizon 
grow globally. This non-local behaviour is encoded in the elliptic nature of Eq. (6), 
providing an example of the non-local behaviour of these dynamical trapping horizons. 
This is possibly disturbing, if considering as the boundary of a physical object. An 
even more curious behaviour of follows from the two following results: 

i) Foliation uniqueness [19]: the foliation by MOTSs of a dynamical FOTH is unique. 

ii) Existence of DHs [12, 13]: Given a (stable) marginally outer trapped surface y§ in 
a Cauchy hypersurface Eo, to each 3+1 spacetime foliation (E ? ) ?g k containing Eo it 
corresponds a unique adapted dynamical FOTHs that contains y§ and is sliced 
by marginally outer trapped surfaces {y} such that y t cE,. 



The first result provides a sort of rigidity for DHs, very useful in our context since it 
determines the evolution vector h a uniquely up to a time reparametrization. Regarding 
the second result, this is a crucial benchmark in the treatment of quasi-local horizons 
in an Initial Value Problem approach. The main point we want to underline here is that 
the combination of results i) and ii) above leads to the non-uniqueness of DHs: the 
combination of results on evolution existence and foliation uniqueness for DHs implies 
the generic non-uniqueness in the evolution of a FOTH from an initial MOTS. To see 
this it suffices to consider the evolution of an initial S^q E Zq into DHs Jt\ and ^2 
compatible with (generic) 3+1 foliations and {E ?2 }, assume Jt\ = J%2 an d then 

reason by contradiction with the result on the uniqueness of the foliation (see e.g. [9]). 
Such non-uniqueness in the evolution is a surprising behaviour if we consider Jf? as the 
boundary of a physical object. Actually such behaviour is a characteristic signature of 
gauge dynamics. At this point it is worth to remark that the amount of gauge freedom in 
the stationary and the dynamical cases remains the same: whereas i) in equilibrium the 
geometric hyperfurface Jf? is unique, but its foliation by MOTS is not unique and can 
be parametrized by a free function / on ,Y rescaling the null normal (£ a — > f£ a ), ii) in 
the dynamical case the geometric hypersurface is non-unique but the foliation by MOTS 
is rigid. In this latter case, the gauge freedom is given by the choice 3+1 foliation, i.e. by 
the lapse function N evaluated on y. Therefore, both equilibrium and dynamical cases 
present the gauge freedom (the choice of a function on y), but dressed differently. 



A pragmatical view on FOTHs 

From the discussion above we must retain: i) FOTHs are hypersurfaces tracking the 
BH region, ii) they are well-adapted to the 3+1 Initial Value Problem approach, and iii) 
DHs incorporate a sort of rigidity fixing the evolution vector up to time reparametriza- 
tion. These are remarkable geometric properties. If considered as physical surfaces, 
FOTHs present curious properties: i) non-uniqueness, ii) superluminal behaviour (when 
dynamical) and iii) global behaviour. The last point has been addressed in detail in Ref. 
[20], where clairvoyance properties of trapped surfaces are discussed, in particular their 
capacity to enter into flat regions, one of the reasons to abandon event horizons. 

Although dynamical trapping horizons have proved very useful to gain crucial insight 
into physical aspects such as BH thermodynamics in dynamical contexts [15, 11], for 
the reasons discussed above we adopt here a perspective underlining the role of FOTHs 
as purely geometric probes to be employed the analysis of dynamical spacetimes. When 
considering an inner test screen J^f( nn in the context of the cross-correlation approach, 
we look for a hypersurface that: i) should be a footprint of the BH presence, in partic- 
ular providing a probe into their spacetime geometry, ii) should be suitable for the Ini- 
tial Value Problem approach, providing in this setting a preferred geometrically defined 
structure with some sort of rigidity, and iii) should intrinsically incorporate an evolution 
concept, tracking the evolution of the BH properties. In this context, we adopt a prag- 
matical approach on DHs, as hypersurfaces of remarkable geometric properties in BH 
spacetimes, providing preferred geometric probes into the BH spacetime geometry. 



FIGURE 2. Illustration of the use of the Initial Value Problem approach in the cross-correlation scheme. 
A DH is used as inner screen ^ nn , whereas the large spheres world-tube 38 (or J^ + ) is used for ,3^ oul . 

Initial Value Problem and cross-correlation approaches 

The adoption of an Initial Value Problem approach to the construction of the space- 
time has an impact on the cross-correlation scheme sketched above. On the one hand, as 
already discussed, the event horizon $ is not available during the evolution. In addition, 
in a Cauchy evolution of an asymptotically flat spacetime one does not constructs null 
infinity J^ + , but rather the slices E f reach spatial infinity r (as an alternative, one could 
solve a hyperboloidal Initial Boundary Value Problem, rather than a Cauchy one, to con- 
struct J f+ ). On the other hand, advanced and retarded times u and v are not natural, and 
one rather employ a 3+1 time function t associated with the slicing {E f }. 

Therefore, in our 3+1 treatment of the cross-correlation methodology we shall employ 
(cf. Fig. 2): i) dynamical trapping horizons as canonical inner screens associated with 
the 3+1 slicing with lapse function N, ii) a timelike worldtube £3 at large distances as 
outer screen (or ^ + in a hyperboloidal slicing), and iii) a 3+1 spacetime slicing time 
function t that automatically implements a (gauge) mapping between u and v. 

In the context of the cross-correlation of quantities h- mn (t) and h out (t) as time series, 
we note that generically the 3+1 slices Z t intersect multiply the dynamical trapping hori- 
zon (cf. Fig. 2). This is the underlying reason for the jumps occurring in the evolution 
of apparent horizons, generic in 3+1 BH evolutions. From the 3+1 perspective both an 
external and an internal Jtf/ m apparent horizon are present. In numerical simula- 
tions it is standard to neglect the internal horizon as irrelevant. The splitting of the single 
spacetime hypersurface ffl 3 in two parts, ffl — J£f nt U J^ xt is an artifact of the 3+1 de- 
scription and simply reflects that the coordinate t is not a good label for M '. If we are 
interested in cross-correlations only after the moment t c of first appearance of the appar- 
ent horizon, then using the external part J^ x t is enough. But if integrating fluxes in time 
is relevant in our problem, then the internal J^ nt must be kept into the picture in order 
to account for the whole history of the flow into the BH singularity (cf. Fig. 3). 



FIGURE 3. Internal and external apparent horizons in a 3+1 BH evolution. After the moment t c of first 
appearance, slices E, intersect multiply the DH Jff, splitting it into Jf eKt and Jtfnt (dashed line) parts. 

APPLICATION TO BLACK HOLE RECOIL DYNAMICS 

We apply now the previous ideas to the study of the recoil dynamics of the BH resulting 
from the asymmetric merger of two BHs. In an asymmetric binary merger the emission 
of linear momentum through gravitational waves is not isotropic, so that the final rem- 
nant must recoil in order to preserve the total linear momentum. This is astrophysically 
relevant in the context of the merger of supermassive BHs in galaxy encounters. 

We aim at gaining insight into the dynamics controlling the recoil (kick) velocity 
Vk and, in particular, into the systematics of a late-time deceleration referred to as the 
anti-kick in the literature. On the one hand, the very presence of such final deceleration 
is a direct consequence of the fact that the kick velocity is obtained upon integration 
of a decaying oscillating quantity, namely the flux of linear Bondi momentum at null 
infinity: there is no reason to expect that the maximum and the asymptotic value of the 
time integral of a decaying oscillating signal coincide [21]. The relevance of this problem 
lays, rather than on the presence itself of the anti-kick, on the capability to estimate a 
priori its magnitude from an understanding of the underlying dynamics. In particular, 
a relative large antikick reduces the resulting kick velocity, whereas large final recoil 
velocities are associated with small relative antikicks. In this context, the magnitude 
of the antikick is controlled by the ratio between the oscillation timescale, T, and the 
decaying timescale, T, leading to the crucial notion of slowness parameter [22] 

T 

P=-, (7) 

so that for P <C 1 the rapid oscillations in the signal induce cancellations in the time 
integral and lead to large antikicks, whereas for P rs 1 the antikick is small. 

Understanding the involved oscillating and decaying dynamics is therefore crucial in 
this problem. Ref. [23] paved the way to get insight into the responsible gravitational 
dynamics, in terms of the analysis of the evolution (dissipation) of the quasi-local 



horizon geometry. The cross -correlation scheme [8, 9, 10] permits to develop a more 
systematic analysis along those lines. Before proceeding further, and since it will be 
relevant later, we note that quasi-local BH horizon tools have already been applied to this 
problem. In particular, Ref. [24] proposes an expression for the BH linear momentum 

m = -Lf {K ab -Ky ab )^s b dA, (8) 
87T Jy t 

by applying the ADM prescription at spatial infinity f to the apparent horizon. 



Cross-correlations in BH recoil dynamics 

In the application of the cross-correlation approach to BH recoil dynamics, we take as 
quantity h out the flux of Bondi linear momentum at (actually at the approximation 
given by large spheres worldtube g$) along a spatial direction t, a . That is 
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dt 



(9) 



where ^V{t) is the so-called news function. At the inner DH screen M', we lack a 
geometric analogue to dP B [^]/dt, in particular no news function formalism is available. 
However, Eq. (9) suggests a natural heuristic candidate K[t,\ for the quantity h mYl 

K[£](t) = ~ I J^\t) 2 dA , A>it) = r^(t')dt' , (10) 

IbTt Jy, J to 

where the Weyl scalar plays the analogous role at ^f nn that ¥4 plays at (the 
superindex N refers to a choice of null tetrad adapted to the 3+1 foliation with lapse 
function N [9]). However, the function ,jVy{t) does not satisfy a key requirement on the 
news function, namely its characterization purely in terms of the geometry of a section 
Sf t of the horizon. In other words, the flux defined by the square of the news must be 
an instantaneous quantity defined by quantities crossing the horizon at a given time. 
This issue can be addressed by correcting the integrand in ,jV^>{t) with terms completing 
^oit') to a total time derivative. From the tidal equation for the evolution of the shear 

ojfi along Jf, with leading order given by *¥q, we propose the more geometric quantity 

The notation dP { ' 0) [%]/dt is meant to underline that this quantity is well defined at ,Y t , 

without implying the existence of a well-defined conserved quantity P 

The idea now would be to cross-correlate quantities (dP [£]/dt)(t) at Jf? and 
(dP^lQ/dt)^) at J r+ . This requires however the determination of the evolution vector 
h a , which involves the resolution at each time step of the elliptic Eq. (6) for C. Although 



no conceptual issues are involved in this, we adopt in a first stage a technically simpler 

strategy where for h mn , rather than dP [g]/dt, we use an effective curvature vector 
constructed in terms of the intrinsic geometry of =5^. In order to justify this, we write 

8 h 2 R = -0 W 2 R + 2 2 D a2 D b o^ - 2 A0 W , (12) 

for the evolution of the Ricci scalar 2 R of the induced metric q ab on 5? t . To fully control 

its evolution, we need to track the evolution of q ab , 9^ and of b . In order to get insight 
into the involved quantities, we make explicit the system for a null horizon (the actual 
spacelike case has the same structure, but involving corrections on the function C). Then 

8?R = -6^ 2 R + 2 2 D a2 D b a^ - 2 AflW 
Stqab = 2o$ + eW qab 

See® = -\{e^) 2 -a^a^ ab -^T ab ti b (13) 

= ^a^ cd qab - q \ q d b C lcf d l i f , 

where C a bc d denotes the Weyl tensor. Once initial data are given, the whole system is 
driven by the external forces given by q c a q d b Ci c fdl l £f and T ab £ a £ b . Focusing here on 
the vacuum case (see [9] for a more general discussion) we note that q c a q d b Ci c fd£ l £f = 
^omafnt, + mom a m b , where m a is a complex null vector tangent to 5? t - On the one hand, 
the evolution of the whole system (13) is determined by the at the horizon, which 

justifies the understanding of ,jV ab in Eq. (11) as a kind of news-like function. On 
the other hand, we note that the evolution of 2 R is completely driven by the rest of the 
system, without backreacting on it. This last point is crucial, since the evolution of 2 R 
then captures in an effective way and in a single function the evolution of the whole 
system, in particular the evolution of the shear (but also other degrees of freedom if 
matter is present). This leads to the introduction of the effective curvature vector [8] 

K^](t) = - -L / {g Si )\J>{t)\ 2 dA , J>{t)= f 2 R(t')dt' + ^ , (14) 
ion Jy, J t c 

as an effective estimator of the evolution of the horizon geometry (here ,jV tc is an 
initial value function; see [8] for details on its fixing). In a first stage of our quantitative 
analysis, we use K [£,] as the quantity h irm to be cross-correlated with dP B [%]/dt. 

In order to test these tools, we have considered the head-on collision of non-spinning 
black holes with mass ratio q = 1/2 and have constructed numerically the associated 
dynamical spacetime (cf. details in [8]). We extract the timeseries corresponding to 
K efl [%](t), once the common apparent horizon has formed and (dP B [^]/dt)(t) at (an 
approximation to) null infinity during the whole evolution. From a qualitative perspec- 
tive, both timeseries show a good agreement, from the moment of first appearance of the 
common apparent horizon. The qualitative agreement is preserved in time (see Fig. 4). 

For a quantitative comparison we use the cross -correlation of h\ (t) and 



h\(t + x)h 2 (t)dt , 
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FIGURE 4. Effective curvature A^ ff (f ) at the horizon from the moment of formation of a common hori- 
zon (red dashed curve) and flux of Bondi linear momentum (dP B [§] jdt) (t) evaluated at an approximation 
of J* + (blue dotted solid curves, before and after the appearance of a common horizon, respectively). 
Panel (a) shows the good qualitative agreement between both quantities after the merger. Panels b) and (c) 
compare the same quantities for latter times, showing the persistance in time of the good agreement. 



which encodes a quantitative comparison between the two timeseries as a function of 
the time-shift % {lag) between them. In particular, the matching between the two signals 
can be expressed in terms of the number 

Jt(h u h 2 ) = max I Corr(h u h 2 )(r) (J6) 



x 



[Corr(/*i , hi) (0) ■ Corr(/* 2 , h 2 ) (0)] : 



This number is confined between and 1 (with 1 indicating perfect correlation, and 
no correlation at all) and provides the maximum matching between the timeseries. The 
calculation of the correlations in our scheme requires a careful treatment of what can 
be referred to as a time stretch issue, in order to deal with the freedom in the choice of 
spacetime foliation that determines the gauge mapping between retarded and advanced 
times u and v (cf. [8] for details). Once this is taken into account, the calculation of 
the correlation parameter in our problem gives typically values > 0.9, independently 
of the width of possible time windows applied to the signals prior to the calculation 
of the correlations [8]. In order to assess the potential bias in the calculation of the 
correlations due to the time decay of the signal, we model the signals by exponentially 
decaying functions, h mn (t) = e~ Kinnt h* nn (t) and h oul (t) = <? _,w /Zo Ut (f), and perform the 
correlation analysis in the time series h* nn (t) and h* ut (t) (cf. left panel in Fig. 5). Again, 
the correlation number j$ « 0.9. More interestingly, Fourier transforming the signals 
to get the power spectrum, we find that only two frequencies enter the dynamics in this 
head-on collision, Q\ nn and Qf n at J4? and Q? ut and Q^ ut at S+ (cf. right panel in Fig. 5). 
From the very good approximation [8] given by Kf f ~ ^¥iJYi and (dPf jdt) ~ JV 2 <M>, 
where jYi and jty are the (m = 0) spherical harmonic components of jV and jY , 
respectively. Remarkably, using a sinusoidal Ansatz for jV^ and JVi, we can reconstruct 
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FIGURE 5. The left panel shows the signals Kf(t) and (dP B [%]/dt)(t) shown in Fig. 4, once the 
exponential decay has been eliminated. The right panel presents the power spectrum of the corresponding 
signals showing clearly the presence of two dominating frequencies il 1 "" at JSf? and i2° U 2 at J? + . 



from at ^ an d ^2 at ^ > tne corresponding 1 = 2 and £ = 3 modes 

oinn r»inn . oinn 1 oinn C»oul oout Ciout 1 r\out 

* L l=2 — 2 ' f = 3 ~~ 2 ' ~~ 2 ' ~~ 2 

(17) 

leading to an extremely good agreement with the real part of the BH quasi-normal mode 
frequencies (of =2 and 0^3- On the other hand, the decay inverse time scales Ki nn and 
Ko Ut can be retrieved from the addition of the imaginary part of the £ = 2 and £ = 3 quasi- 
normal modes: K^ecay = ®} = 2 + a) i=3- This matching of the frequencies at Jf? and J + 
with the quasi-normal modes is shown in Table 1 . On the one hand, this identification 
of the role played by the quasi-normal modes is consistent with the simple dynamics 
of the gravitational field in vacuum. On the other hand, it impacts directly our specific 
recoil dynamics problem since characteristic decay and oscillation timescales can be 
constructed from the imaginary and real parts, respectively, of the quasi-normal modes 

2% 2% 

T= „„ , „, , (18) 



leading to the slowness parameter introduce in Eq. (7) 

T co} =2 + co} =3 



T (Of, + co; 



(19) 



This expression has a predictive power to estimate the recoil of the final BH remnant, 
from an initial configuration of binary BHs. Indeed, using analytic estimations of the 
final BH M mass and spin parameter a corresponding to given initial configurations 
[25], one can calculate the associated Kerr quasi- normal modes and construct P in 
(19). Further insight into this slowness parameter P is gained from the dynamics of 
the geometry of the dynamical horizon JF. As shown in Eq. (12), the expansion flw 



TABLE 1. The first row shows the comparison between the £ = 2,3 oscillation fre- 
quencies at J4? (i.e. 3 ) and at (i.e. £i^ 2 3 ) with the real part of quasi-normal 
modes of a Schwarzschild BH, (of =2 3 . The second row shows the comparison between 
the decay exponents K"j nn and K^, ut , respectively at and J? + , with the decay coefficient 
calculated from the imaginary part of the quasi-normal modes as JQjecay — a> e=2 + (0 t=3- 
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0.60 ±0.04 


0.59 ±0.04 0.59944 
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0.181 ±0.006 


0.179 ±0.005 


0.18166 







controls the dynamical decay of the intrinsic geometry. This dissipative role of 6^ is 
further supported by its interpretation in the membrane paradigm [26, 27, 28, 29, 16, 30, 

(h) 

31, 32] (see also discussion in [9]) as associated with bulk viscosity terms, whereas o\ 
is related to the shear viscosity. The latter is responsible for the (shape) oscillations in the 
intrinsic geometry (cf. last equation in (13), for its relation to propagating gravitational 
degrees of freedom encoded in the Weyl tensor). Given the physical dimensions [9^] = 
[ctW] = [Length]" 1 , one could introduce decay and oscillation inverse timescales by 
averaging 9^ and <j( h \ respectively, over the horizon section y t . In order to make 
more precise this heuristic approach, we consider the evolution equation for 9^ [31] 

(fc+ew)©**) = -^flW+^ffW (20) 

+ (9^Y_ + 2 Da( 2 DaC _ 2C£l\P) + SnT ab T a h b - 9^8 h C (21) 



Denoting by the unit vector in the instantaneous spatial direction of motion of the 
BH, decay and oscillation timescales can be introduced as 

1 1 



= ^j^ Si )(K^9^)dA 



x{t) 2 A Jy, 

1 1 r . ( .(h) Mab 



so that the instantaneous slowness parameter 



K> W ) dA > (22) 



n _ T(t) I §y&*i)(* h) e {h) )dA ( 

satisfies P(t) « 1 near equilibrium, when derivative and higher-order terms in (20) can be 
neglected, only surviving the terms in the first line of the right-hand side, precisely those 
used to define T(t) and x[t) (note that they lead to the Hartle-Hawking area evolution 
equation). Note that P ps 1 is consistent with the absence of antikick near-equilibrium. 



Contact with quasi-local BH linear momentum 



Quasi-local notions of linear momentum has been applied to the study of BH recoil 
dynamics in [24, 33]. From this perspective, the news-like function jY^ in Eq. (11) 
can be used to introduce a heuristic flux of Bondi-like momentum at the horizon M '. 
Consider the (timelike) unit normal x a = t= = -k^(£ a + Ck a ) = -j=(bn a +Ns a ) to 

Jft? (here, n a is the timelike unit normal to a 3-slice L t and s a is the spacelike unit normal 
to y t in L t ). Consider also a generic 4- vector r\ a . We can introduce a flux of Bondi-like 
4-momentum as [9] (note the natural use of an advanced time parameter v) 

T - -W™ - -^jm, (< ^>> . ( 24) 

Considering an Eulerian observer n a , the associated flux of energy would be 

dP s dPM = lfb ( a (%V.)°>) dA , (25) 
dv v ' dv 167tJyy/2C\ ab ' 

with = a/1 +N 2 /2C, whereas the flux of Bondi-like linear momentum is given by 

dv 16k J.y v v^C V ^ ' \ ab ) 

These expressions are closely related to those proposed for DHs [17, 18]. The integration 
in time v along the horizon ffl provides a heuristic prescription for the linear momentum, 
a sort of Bondi-like counterpart of the heuristic prescription in Eq. (8) based on the ADM 
momentum. More generally, for a quantity Q(v) with F(v) flux through y, we write 

Q(v) = Q(v )+ fF(v')dv' . (27) 

JVQ 

Note that, in view of the 3+1 description (cf. Fig 3) such an expression can be split into 

Q(t) = Qo+ f (F) mt (t')dt'+ f (F) ext (t')dt' + Res(t) 

Jt c Jt c 

where (F) int is the flux through the internal horizon, (F) ext is the flux through the ex- 
ternal horizon and a residual Res(?) = F m (t')dt' must be included (a more complete 
expression taking into account changes in the metric type of is discussed in [9]). This 
expression makes explicit the relevance of tracking the internal horizon when addressing 
the integration in time of physical fluxes across the dynamical horizon Jf. 



CONCLUSIONS AND PERSPECTIVES 



We have outlined some basic elements of a cross-correlation approach to the analysis 
near-horizon spacetime dynamics. In particular, we have identified DHs as hypersurfaces 



providing inner canonical screens in a 3+1 Initial Value Problem approach to the space- 
time construction, where appropriate geometric quantities can be defined to probe and 
monitor bulk dynamics, namely through correlation with quantities at an outer screen. 
We have applied this scheme to the study of BH recoil dynamics. First, we have in- 

troduced a heuristic horizon news-like function to build a quantity dP [g]/ dv tracking 
quasi-locally (on Jt?) the qualitative and to a good extent also the quantitative features of 
the recoil dynamics at In particular, the analysis of its correlations with the flux of 
Bondi momentum at in numerically constructed spacetimes corresponding to binary 
BH head-on collisions confirms the simple character of vacuum spacetime dynamics in 
this setting. Second, this latter remark has led to the proposal of a prescription for a 
slowness parameter P, as constructed from the complex BH quasi-normal modes. More 
generally, inspired by the BH horizon viscous fluid analogy in the membrane paradigm 
and further supported by the horizon geometry dynamics, we have proposed geomet- 
ric decay and oscillation timescales leading to a more general characterization of the 
slowness parameter P, with a well-defined instantaneous meaning. Third, we have made 
contact with the heuristic attempts for estimating the BH 4-momentum, with the pro- 
posal of a quasi-local Bondi-like expression. In this context, we have emphasized the 
relevance of the internal 3+1 horizon when considering flux integrations along the BH 
horizon Jtf. An open problem in this sense is to assess the capability of a DH to dress 
the BH singularity, this involving the understanding of early and late DH asymptotics. 

Despite the encouraging prospects, there are important caveats in the sketched cross- 
correlation approach. First, a sound proper formalism is missing. In particular, an "in- 
verse scattering" picture must still be systematically developed. In addition, the gauge 
issue briefly referred to as the "time stretch issue" must be addressed in generic situ- 
ations. More importantly, it is not clear how to assess the conditions under which the 
comparison/cross-correlation of quantities at outer and inner screens is actually legit- 
imate. A possible approach to these issues (cf. [10]) would consist in considering the 
cross-correlation of test-fields evolving on dynamical spacetimes, without backreacting 
on them, so that the field evolution faithfully tracks specific relevant aspects in the ge- 
ometry of the spacetime dynamics (see also [34]). On the one hand, this would remove 
the ambiguities in the choice of quantities h mrx and h out to be correlated at ^f nn and 
J^out- On the other hand, it would permit to extend the scheme to the direct analysis of 
bulk quantities (this would be related to the techniques in [35, 36] for tracking spacetime 
dynamics). In particular, the strategy of analyzing the spacetime geometry in terms of 
correlations of appropriate test-fields, that could be paraphrased as pouring sand on a 
transparent surface, can benefit from the use of tools and concepts developed in statis- 
tical approaches to field theory in curved spacetimes, e.g. [37, 38, 39, 40, 41, 42, 43]. 
This attempt leads to the following declaration of intentions as a perspective for future: 
to develop a strategy for spacetime analysis aiming at a functional and coarse-grained 
description of the spacetime geometry, by importing functional tools for the analysis of 
condensed matter and quantum/statistical field theory systems ( in curved backgrounds). 
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